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In a recent paper Dajka, et.al., [J. Phys. A 40, F879 (2007)] predicted that some 
composite systems can be entangled forever even if coupled with a thermal bath. We an- 
alyze the transient entanglement of a single-Cooper-pair box biased by a classical voltage 
and irradiated by a quantized field and find the unusual feature that the phase-damped 
cavity can lead to a long-lived entanglement. The results show an asymptotic value of the 
idempotency defect (concurrence) which embodies coherence loss (entanglement survival), 
independent of the interaction development by dependent critically on environment. 



1 Introduction 

Josephson junctions are being investigated as a possible route to scalable quantum com- 
puters [1]-[TT]. The present lack of a current standard based on quantum devices has 
inspired several attempts to manipulate single electrons, where the rate of particle trans- 
fer is controlled by an external frequency. One of the physical realizations of a solid-state 
qubit is provided by a Cooper pair box which is a small superconducting island con- 
nected to a large superconducting electrode, a reservoir, through a Josephson junction 
[T2] . Also realizations of superconducting charge qubits are a promising technology for 
the realization of quantum computation on a large scale p!3]-|17]. 

In this context, a solid-state system is highly desirable because of its compactness, 
scalability and compatibility with existing semiconductor technology. Even though a 
Cooper pair box can contains millions of electrons at any one time, the box exhibits only 
two quantum charge states, depending upon whether or not a Cooper pair of electrons 
has recently tunneled into the box and various superconducting nanocircuits have been 
proposed as quantum bits (qubits) for a quantum computer [9l|10]. By gating the Cooper 
pairs into the box with an appropriate pulse width, previous research has shown that 
a coherent superposition of the two states can enable quantum computations. In archi- 
tectures based on Josephson junctions coupled to resonators, the resonators store single 
qubit states, transfer states from one Josephson junction to another, entangle two or more 
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Josephson junctions, and mediate two-qubit quantum logic. In effect, the resonators are 
the quantum computational analog of the classical memory and bus elements. 

The present work is motivated by conjectures and statements presented in a recent 
fast track paper [18J and experimental results on Josephson junction and normal metal 
flux qubits coupled to the environment [6]. We obtain a long-lived entanglement using 
a superconducting charge qubit. More precisely, we endeavor to show the important 
property of entanglement via idempotency defect of a single Cooper pair box, due to the 
presence of a phase-damped cavity. Despite the complexity of the problem, we obtain 
a quite simple exact solution of the master equation that is valid for arbitrary values of 
the phase damping. In the framework of the exact solution of the master equation, we 
determine the coherence loss and the degree of entanglement. We perform a systematic 
analysis in order to reach an understanding of the Cooper pair dynamics in the presence of 
the decoherence. Physically, the effect of phase damping may be understood to be analogy 
of the T2 spin depolarization effects observed in nuclear magnetic resonance spectroscopy 
(for detailed physical motivation see Ref. [7J). Besides phase-damping-model importance 
in the description of different physical situations, it is very instructive since it allows for 
obtaining analytical treatments for different entanglement measures of some classes of 
states [8]. Some theoretical discussions and analysis of special cases of the problem at 
hand were given in Refs. [H dSl [20] and experimental results were predicted in Ref. [H] . 

The organization of this paper is as follows: in section 2 we introduce the model and 
give the exact solution of the master equation. In section 3 we employ the analytical 
results obtained in section 2 to discuss the idempotency defect and entanglement for 
different values of the phase-damped cavity. Finally, we summarize the results in section 
4. 

2 The model 

Several schemes have been proposed for implementing quantum computer hardware in 
solid state quantum electronics [21\. These schemes use electric charge, magnetic flux, 
superconducting phase, electron spin, or nuclear spin as the information bearing degree 
of freedom. 

We start our analysis by presenting a brief discussion and a few physical principles of 
the Cooper pair box system. We consider a superconducting box with a low-capacitance 
Josephson junction with the capacitance Cj and Josephson energy Ej, biased by a clas- 
sical voltage source Vg through a gate capacitance Cg and placed inside a single-mode 
microwave cavity. In particular, the schematic picture of this single-qubit structure may 
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Figure 1: Schematic picture of the Cooper-pair box which is driven by an apphed voltage 
Vg through the gate capacitance Cg. Black bars denote Cooper pair box. The two Joseph- 
son junctions have capacitance Cj and Josephson energy Ej. The driving microwave field 
is generated using the electrical voltage acting on the charge qubit via the gate capacitance 



be modeled as shown in figure 1. The total Hamiltonian of the system can be written as 
[22] 

H = \{Q- CgVg - CgVf{Cg + C j)'^ ~ Ej COS + nwi^^tfj + 

where Q = 2Ne is the charge on the island (e is the electron charge and N is the number 
of Cooper-pairs) and is the phase difference across the junction. The radiation field 
is to produce an alternating electric field of the same frequency across the junction, and 

V is the effective voltage difference produced by the microwave across the junction. We 
assume that the dimension of the device is much smaller than the wavelength of the 
applied quantized microwave (which is a realistic assumption), so the spatial variation 
in the electric field is negligible. We also assume that the field is linearly polarized, and 
is taken perpendicular to the plane of electrodes, then V can be written as [IHl [23] 

V = ifkj (^ip — /{2Cf), where i/'^ and i/j are the creation and annihilation operators 
of the microwave field with frequency u. We denote the capacitance parameter by Cp, 
which depends on the thickness of the junction, the relative dielectric constant of the thin 
insulating barrier, and the dimension of the cavity. 

We consider the case where the charging energy with scale = {Cg + Cj) , which 
dominates over the Josephson coupling energy Ej and concentrates on the value Vg = 
e/Cg, so that only the low-energy charge states = and = 1 are relevant. In this 
case the Hamiltonian, in the basis of the charge states |0) and reduces to a two-state 
form. In a spin- 1/2 language 



H = E,{1 + e-^C'jV) - ^Eja, + 2e~'E,CjVa, + hw (^^V + ^ , (1) 

where a^c and are the Pauli matrices in the pseudo-spin basis. It is to be noted that 
the charge states are not the eigenstates of the Hamiltonian ([Tj), so the Hamiltonian can 
be diagonalized yielding the following two charge states subspace |e) = — |0)) and 

\g) = + 10)). Here we employ these eigenstates to represent the qubit. If we consider 

a weak quantized radiation field and neglect the term containing V"^, the Hamiltonian ([T|) 
can be rewritten in the rotating wave approximations as 

H = M,^^ + \e,.. + { j^I^y^,)., + H.C.] . (2) 

We consider the interaction with an environment to be as the phase-damping type. This 
is a reservoir coupled to the field via the number operator of the indicating field, so that 
there is no energy damping, although there is a phase damping. 

In order to obtain the general solution of the master equation for the density matrix 
under the phase damping of the cavity field at a zero temperature bath, we write 

^ = ^[i^. Pit)] + 7 (2V^Vp(t)^V - ^V^Vp(t) - pW^V^^^) , (3) 



where 7 is the phase-damping constant. An important feature of this quantized system is 
that its steady states, known as dressed states, are entangled. Switching to an interaction- 
picture representation for convenience by defining p'{t) = exp{iHt)p{t) exp{—iHt) , exact 
solution of equation (4) can be obtained in the dressed-states representations [25]. Con- 
sequently, equation ([3]) can be written as 

^^fM = ^e^^* (^2V'VpWV'V - V-V^^^W - p(t)^V^V) e-'""'. (4) 

Next, we write the field operators and ^p'ip in terms of the dressed states basis and get 
the initial state of the system expressed in the product density matrix forms. In the basis 
In, e) and \n, g) states, the field operator ^ can be written as ■)/' = ^ y/n{\n — 1, e){n, e\ + 

n=0 

\n — 1, g){n, g\). Neglecting the oscillating terms of the master equation (jl]) in secular 
approximation, the density matrix in terms of the dressed states becomes 

n=0 

+ I0i^^)(0^^|p ((2^ + 1)' + 1) + I0^^)(0i~^|p W 

00 

+1 {((2^ + l)(2^ + l)(0it-VW+PW0L"^)}, (5) 

n,m=0 

where \(f)^^) are two eigenstates of the Hamiltonian (3) for a lossless cavity, 0^'tn~^ — 
\<P^n^){(p^n^\ + 10m ^) (0m ^ I ^nd = /i„ - The eigenvalues are given by ±/x„, where 



We denote by A = Ej — uj the detuning between the Josephson energy and cavity field 
frequency, {5 = A/2). Based on the preparatory work, now we can find an exact solution 
under certain conditions of the whole system. With this in mind we will assume that the 
initial state is prepared to be a particular coherent state P'^{0) = |a)(a| of the field with 
the Cooper pair box prepared in the state p"^(0) = |e)(e|. The initial state of the system 
can be expressed in the product density matrix form, p(0) = p'^(0)(S) p-^{0)- Consequently, 



the general solution to equation ([5]) my be written explicitly as 



pit) = ^"^^ ( ~Y^] {-^t{n - m)^) 



n,m=0 




X |exp(-i/5i2) (cos (/i„„(t)) + cos (/i'„^(t))) |n, e) (m, e 
"2 exp(-^/3l2) sin |n, e) (m + 1, ^| 

+ - exp(i/5i2) sin {fi,^^{t)) \n + 1, ^) (m, e| + exp(-i/5i2) 
X (cos (/i„„(t)) - cos (^'„„(t))) |n + 1, g) (m + 1, ^1 1 , 



(7) 



where finmit) = Pni^) — Pmit) , p'nmit) = Pnif) + Prn{t) ■ The probability distribution among 
Fock states is Poissonian, 6„ = {n\a), with n = \a\^ and (3 12 — where (3 is the phase 

of the initial state of the field i.e. a = |a|e*^. The decoherence effect on the dynamical 
evolution of the present system can be discussed through the phase-damping constant 7. 

3 Coherence loss and entanglement 

In general, due to decoherence, a pure state is apt to change into a mixed state. However, 
in many cases of quantum information processing, one requires a state with high purity 
and large amount of entanglement. Therefore, it is necessary to consider the mixture of 
the state and its relation with entanglement. 

Here we use the idempotency defect, defined by linear entropy, as a measure of the 
degree of mixture for a state p'^(t), in analogy to what is done for the calculation of the 
entanglement in terms of von Neumann entropy ^B] which has similar behavior. In order 
to analyze what happens to the Cooper pair box, we trace out the field variables from the 
state p{t) and get the reduced density matrix p'^(t) = trfp{t). The idempotency defect as 
a measure of coherence loss can be written as 



where 8l has a zero value for a pure state and 1 for a completely mixed state. 

Supplemental to the analytical solution presented in the above section, here we discuss 
the results obtained numerically and interesting situations occurring for different values 
of the detuning and phase-damped cavity parameters. We consider the experimental 
parameters, described above, which are accessible using the present day technology as, 
Cj ~ u ^ 10^°Hz, Cf ~ KbT « Ej ^ hw « and the initial state 

of the filed as coherent state. In order to analyze the effects resulting from variation in 
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Figure 2: (a,b) Plot of idempotency defect as a function of the dimensionless scaled time 
At (A = ^Je^uj/ihCp) and the detuning parameter A/A and (c) the atomic inversion as 
a function of the scaled time At. In the contour plot, disentanglement is shown in the 
severe shading 



the detuning or phase damped cavity we consider the idempotency defect as a function 
of the scaled time Xt and A/ A (7/ A) shown in figures 2 and 3. We have fixed the mean 
photon number of the coherent field as n = 10. As can be seen from figure 1, Sj;'^^ 
smoothly diminishes with increasing the detuning parameter. For further increasing of 
the detuning the impurity of the state of the Cooper pair box system is rapidly growing 
and Sj:"^^ disappears completely. For the case when we take 7 = A = 0.0, we get almost 
zero values for the idempotency defect only at t = 0, which means that a pure state will 
not be reached at any time except at the initial stage of the interaction time (see figure 1). 
To apprehend the essential features of detuning effects on coherence loss, we presented in 
figure lb the contour plot of the concurrence, where complete separable states are shown 
in the severe shading 

In Figure 2c we show the time evolution of the atomic inversion. Apparently, it is easy 
to observe the existence of collapse and revival of Rabi oscillations of the atomic inversion 
and the first maximum of the idempotency defect is achieved in the collapse time, while 
at one-half of the revival time, the idempotency defect reaches its local minimum. Also, 
it is noticed that in absence of both detuning and phase- damping, a gradual decrease in 
the amplitudes of the Rabi oscillations is shown. 

On the other hand, the decoherence introduces irreversibility into the Junction dy- 
namics and also on the global system. During the repeated periods of maximum and 
minimum idempotency defect, the states of the Junction and field lose and gain coher- 
ence, but given the continuous amplitude decreasing of coherent states, the coherence 
recovered by the Junction is never that which was lost. We may refer here to the work 
given in reference where engineering maximally entangled states has been discussed 
for different systems. Of course, larger the value of 7, the more rapid is this phenomenon 
in the sense of the idempotency defects being close to one (the purity loss of the junc- 
tion state is complete). In particular, for the limiting case of large 7 (7 = O.IA), the 
idempotency defect blows up from zero and rapidly saturates i.e. as time goes on a long- 
lived coherence loss is observed (see figure 2 (top)). Even in a weak-damping cavity, 
the difference between consecutive local maximum and minimum diminishes with time, 
since idempotency defect tends to asymptotic values. Speaking specifically, it arrives at 
a maximum value (about 1) at large values of the phase damping parameter, and then 
remains nearly invariant regardless of the increase of time or 7, while the idempotency 
defect always remains vanishing at At = 0. 

We can gain further physical insight into the dynamical effect of the phase damping 
by considering the general case (mixed state entanglement). To measure the degree of 
entanglement for mixed states of bipartite systems composed by two-level subsystems. 
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Figure 3: Plot of idempotency defect as a function of the scaled time At and the decoher- 
ence parameter 7/A. In the contour plot, a complete mixture is shown in the non-shaded 



one needs to consider a commonly used measure such as the concurrence [28] which has 
been proven to be a reasonable entanglement measure or negativity [29]. Analysis of 
the entanglement decay rates under decoherence for different models of the interaction 
between systems of arbitrary dimensions with the environment has been presented [S]. 
For the density matrix p{t), which represents the state of a bipartite system, concurrence 
is defined as 

C{p) = max{0, ^1-^2-^3- 3^4}, (9) 

where the are the non-negative eigenvalues, in decreasing order (3?i > 3^2 > 3^3 > 
9^4), of the Hermitian matrix T = \f\f^^\f^ and p = i^y ®ay)p* {ay ®Oy). Here, p* 
represents the complex conjugate of the density matrix p when it is expressed in a fixed 
basis and dy represents the Pauli matrix in the same basis. The function C{p) ranges 
from for a separable state to 1 for a maximum entanglement. 




Figure 4: Development of the concurrence C(p) as a function of the scaled time At and 
the decoherence parameter 7/ A. The other parameters are the same as figure 3. 

In figure (jll), we plot the numerically evaluated results for the concurrence C(p), as 
a function of the scaled time At and phase damping parameter 7 in units of A. As 
confirmed in figure (jl]), the asymptotic value of the concurrence is obtained when the phase 
damping is increased. Of course, there are some differences between the concurrence and 
idempotency defect in the amplitudes but the general behavior is the same i.e. comparing 
figures 3a and 4, one can find that concurrence results in qualitative analogy with the 
results of the idempotency defect. This may be thought to arise from the asymptotic limits 



which have been observed in both figures 3 and 4 due to the phase damping. Although the 
entanglement, as witnessed by the concurrence, is lower than maximal possible (about 1), 
it has a fixed value as the phase damping increased (long-lived entanglement). We have 
confirmed the predictions of this phenomenon using a systematic numerical analysis where 
a number of relevant parameters has been varied. However, once the initial state setting of 
the Cooper pair box is considered as p"'(0) = cos^(6')|e)(e| + siia^{9)\g) {g\, {0 < 9 < 7r/2), 
this feature no longer exists and entanglement vanishes in an asymptotic limit. 

Obviously, the above novel phenomena are directly related to the recent results of 
Ref . |18j . With this at hand, one may envision quantum computers using these long-lived 
entangled states for quantum memory and for extended quantum information processing 
[5U] where, superconducting single Cooper pair boxes using superconducting single elec- 
tron transistor fabricated on the same chip as an electrometer has been presented in Ref. 
[31] and the electronic control of a single-qubit achieved in a solid state device has been 
demonstrated [32j. In these works, it has been shown that, the general scalability of such 
a solid state device will be a prerequisite for a practical quantum computer. Also, it has 
been shown that [30J only twice the resources (qubits + elementary quantum gates in 
the decoherence free subspace) are needed to realize up to 4 orders of magnitude more 
operations before the quantum information is lost to the environment. 

It has been predicted only recently that the one-body and two-body responses to a 
noisy environment can follow surprisingly different pathways to complete decoherence 
[33| . The first experimental work and impressive results in this new domain have been 
reported in Ref. [35]. They have devised an elegantly clean way to check and to confirm 
the existence of so-called entanglement sudden death, a two-body disentanglement that 
is novel among known relaxation effects because it has no lifetime in any usual sense, 
that is, entanglement terminates completely after a finite interval, without a smoothly 
diminishing long-time tail |36l [37] . 

4 Conclusion 

In conclusion, we suggest that applying a microwave field to a Cooper-pair box via the 
gate capacitance, a long-lived entanglement can be realized i.e. the Cooper pair composite 
system is entangled forever. In our work we have extended the exactly solvable model 
of a single-Cooper-pair box model by taking into account the decoherence effect on the 
purity loss and entanglement. Decoherence is a very useful concept that has recently 
been widely investigated and has turned out to be very prolific. It is intuitively related 
to the loss of purity of a final state of the quantum system. However, it is demonstrated 



that, it is not correct to think that a quantum system, by increasing the decoherence, will 
suffer an increasing loss of quantum coherence. It is worth stressing in this respect, an 
appropriate choice of the system parameters, specifically, large values of the decoherence 
parameter and initial state setting of the Cooper pair box does give an interesting effect 
to the entanglement process as a long-lived entanglement which may lead to unexpected 
applications. 

We are sure that our ground breaking work on the dynamics of quantum entanglement 

in the Cooper pair box system exactly, will lead both to understand the generic behaviors 
of theses systems by model studies and to add more features to the theoretical models that 
can provide a closer depiction of reality, captured in the near future by higher precision 
experiments. A topic that remains open in almost all decoherence discussions, however, 
is the preservation or destruction of two-body quantum coherence when both bodies are 
small. We are convinced that future experiments exploiting the particular advantages of 
these models will reveal interesting new phenomena and show many surprises. 
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